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Abstract 

Prompt photons are a powerful tool to study heavy ion collisions. Their production rates provide access to the 
initial state parton distribution functions and also provide a means to calibrate the expected energy of jets that are pro- 
duced in the medium. The ATLAS detector measures photons with its hermetic, longitudinally segmented calorimeter, 
which gives excellent spatial and energy resolutions, and detailed information about the shower shape of each mea- 
sured photon. This provides significant rejection against the expected background from the decays of neutral pions in 
jets. Rejection against jet fragmentation products is further enhanced by requiring candidate photons to be isolated. 
First results on the spectra of isolated prompt photons from a dataset with an integrated luminosity of approximately 
0.13 nb~' of lead-lead collisions at -\Z^nn = 2.76 TeV are shown as a function of transverse momentum and centrality. 
The measured spectra are compared to expectations from perturbative QCD calculations. 



1. Introduction 

While di-jet asymmetry distributions at the Large Hadron Collider (LHC) [fl] are consistent with energy loss in 
the hot, dense medium, the detailed physical mechanism of "jet quenching" is still not understood. One of the limiting 
factors in understanding jet quenching is having a proper calibration of the initial jet energy. Replacing one of the 
jets with a penetrating probe, such as a photon or electroweak boson (W or Z), offers the possibility of calibrating the 
energy of the initial jet. This was first proposed by Wang and collaborators in Ref. 0. 

Prompt photons are themselves expected to arise via direct emission and fragmentation from jets. In order to 
remove background events from di-jet processes, as well as the fragmentation photons, an "isolation" criterion is 
typically applied within a cone of a well-defined radius relative to the photon direction. In lower-energy heavy 
ion collisions, the PHENIX experiment performed measurements of direct photon rates in gold-gold collisions at 
y^NN = 200 GeV 0. Tli e CMS experiment performed the first measurement of isolated prompt photon rates in both 
proton-proton as well as lead-lead collisions, at the nucleon-nucleon center-of-mass energy of ysNN — 2.76 TeV |4|. 
They found rates consistent with a scaling with the number of binary collisions. 

This work presents the measurement of the yield of prompt isolated photons in the transverse momentum range 
Pj = 45 - 200 GeV as a function of collision centrality from the ATLAS experiment. The results are described in 
more detail in Ref. 0. 



2. Experimental setup and photon reconstruction 

The ATLAS detector is described in detail in Ref. (6). The ATLAS inner detector is comprised of three major 
subsystems: the pixel detector, the semiconductor detector (SCT) and the transition radiation tracker (TRT), which 
cover full azimuth and pseudorapidity out to \tj\ = 2.5. The ATLAS calorimeter is a large-acceptance, longitudinally- 
segmented sampling calorimeter covering \rj\ < 4.9 with electromagnetic and hadronic sections. In test beam environ- 
ments and in typical proton-proton collisions, the calorimeter is found to have a sampling term for electromagnetic 
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showers of 10-17%/ -\fE\GeV]. The total material in front of the electromagnetic calorimeter ranges from 2.5 to 6 
radiation lengths as a function of pseudorapidity, except the transition region between the barrel and endcap regions, 
in which the material is up to 1 1.5 radiation lengths. The hadronic calorimeter section is located radially just after the 
electromagnetic calorimeter. Within \tj\ < 1.7, it is a sampling calorimeter of steel and scintillator tiles, with a depth 
of 7.4 hadronic interaction lengths. In the endcap region it is copper and liquid argon with a depth of 9 interaction 
lengths. 

To reconstruct photons in the context of a heavy ion collision, the large background from the underlying event 
(UE) is subtracted from each event. This is performed during the heavy ion jet reconstruction, which precedes the 
photon reconstruction, as explained in detail in Ref. [7 1 . This procedure provides a new set of "subtracted" cells, from 
which the mean underlying event, as well as the large-scale modulation from elliptic flow, has been removed. Acting 
upon the subtracted cells, the ATLAS photon reconstruction [8 | is seeded by clusters of at least 2.5 GeV found using 
a sliding window algorithm applied to the second sampling layer of the electromagnetic calorimeter, which typically 
absorbs over 50% of the deposited photon energy. The energy measurement is made using all three layers of the 
electromagnetic calorimeter and the presampler, with a size in the barrel region of Arj x A<p = 0.075 x 0.0125. In the 
dense environment of the heavy ion collision, the photon conversion recovery procedure is not performed, due to the 
overwhelming number of combinatoric pairs in more central collisions. 

To reject clusters arising from hadronic fragments of jets, particularly neutral mesons, the calorimeter is also used 
to measure an isolation energy for each photon candidate, £r(^iso)- The isolation energy is the sum of transverse en- 
ergies in calorimeter cells (including hadronic and electromagnetic sections) in a cone defined by R — tJAtj 2 + A<f> 2 < 
Ri s0 around the photon axis, excluding central core of cells in a region with a size corresponding to 5x7 second layer 
cells. 

3. Data and simulation samples 



The data sample analyzed here is from the 2011 LHC heavy ion run, colliding lead nuclei at -^snn = 2.76 
TeV. After a trigger requiring a 16 GeV energy deposition in the electromagnetic calorimeter, which is sensitive in 
particular to both electrons and photons, events were selected which contained a reconstructed photon or electron 
candidate with a cluster transverse energy of at least 40 GeV. Events were then further analyzed if they satisfied a 
set of quality cuts: the event had to be taken during a period when the detector was found to be working properly, 
leaving an integrated luminoisty of approximately L mt = 0.13 nb~ 1 for this analysis, after the exclusion of several runs . 
Two sets of Minimum Bias Trigger Scintillators covering 2.09 < \rj\ < 3.84 are required to have a well reconstructed 
time signal, and a relative time between the two counters of less than 5 ns. Finally, a good vertex is required to be 
reconstructed in the ATLAS detector, to reject background events from e.g. cosmic rays. In this measurement, yields 
are presented per minimum bias collision, which are estimated using the total integrated luminosity, measured using 
the ATLAS luminosity detectors ||9) . 

The centrality of each heavy ion collision is determined using the sum of the transverse energy in all cells in the 
forward calorimeter (3.1 < \rj\ < 4.9), at the electromagnetic scale. The minimum bias trigger and event selection 
criteria were studied in detail in the 2010 data sample [13] and it was found that 98 ± 2% of the total inelastic cross 
section is sampled by the trigger and event selection requirements. For this analysis, the data have been divided into 
four centrality intervals, covering the 0-10%, 10-20%, 20-40% and 40-80% most central events. In this convention, 
the 0-10% interval has the highest multiplicities, and the 40-80% the lowest. 

For the extraction of photon performance parameters (efficiencies, photon energy scale, isolation properties), a 
set of 450,000 photon+jet events generated using the ATLAS MC11 tune of PYTHIA 6.4 [ 10 1 at yfs = 2.76 TeV, 
is overlaid on minimum-bias HIJING [11] events, which are referred-to as the "PYTHIA+HIJING" sample. The 
generated events are fully simulated using GEANT4 [12] and digitized to produce simulated raw data files, that are 
reconstructed and analyzed exactly as is done for experimental data. 

4. Photon selection and reconstruction performance 

The fine-grained, longitudinally segmented calorimeter utilized in ATLAS allows detailed characterization of the 
shape of each photon shower, providing tools to reject jets and hadrons, while maintaining high efficiency for the 
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Figure 1 : Distributions of photon isolation energy in a i? lso = 0.3 cone for the most peripheral and central centrality bins in data (black points) and 
for MC (yellow histogram), normalized for negative Ej(Ri so = 0.3) values [5]. The differences at large values of £r(Siso = 0.3) can be attributed 
to the presence of jet contamination in the data, which is not present in the MC sample. 



photons themselves. In this analysis, nine shower-shape variables are used, all of which have been used extensively 
in previous ATLAS measurements, particularly the measurement of prompt photon spectra as a function of pseudora- 
pidity l!T4l[T5ll . The cuts used in this analysis are "HI tight" cuts, defined in Ref. J5) as a minimal set of changes to 
the standard set of tight cuts used for unconverted photons in the proton-proton analysis. Photons are also restricted 
to be in the pseudorapidity interval I77I < 1.3. 

Relative to the energy calibration determined in pp collisions [8 1, the residual correction needed to compensate 
for the mixture of converted and unconverted photons is found to be generally less than a percent for much of the 
measured p T range. It was found in the PYTHIA+HIJING sample, that there are no significant differences in the 
residual corrections between peripheral and central events. 

Figure [T] shows the distributions of the isolation energy Ej(R\ S0 = 0.3) for two centrality bins, compared with 
simulated distributions, normalized in the region Ej(Ri so = 0.3) < 0. Both data and MC distributions grow noticeably 
wider with increasing centrality, and agree well in the region of negative isolation energy. The observed enhancement 
on the Ej(R[ so = 0.3) > side of these distributions is expected from di-jet background events. 

5. Measurement of photon yields 

Di-jet backgrounds are subtracted using the so-called "double sideband method" |[T4l [131 ) where photons are 
binned along two axes: the isolation energy within a chosen isolation cone size, and the photon selection criterion 
("tight" or "non tight"). Non-tight photons fail at least one of the more stringent cuts, and have an enhanced di-jet 
contribution. The double sideband approach uses the ratio of counts for non-tight photons which satisfy or fail the 
isolation requirement to extrapolate the measured number of tight non-isolated photons into the signal region. The 
purity extracted using this technique is shown in the left panel of Figure|2](as 1 -Purity) for four centrality intervals. 

The final conversion of the measured yield into a yield per event requires two more factors: the total number 
of events in each centrality sample, and a reconstruction efficiency including all known effects. The efficiencies 
are defined for "HI tight", isolated photons relative to all PYTHIA photons with an isolation energy in a cone of 
R\ s0 = 0.3 around the photon direction of less than 6 GeV. The needed efficiency corrections are categorized into three 
broad classes: reconstruction efficiency, identification efficiency, and isolation efficiency. These are defined in such a 
way that the "total efficiency" e tot is simply the product of these three factors. The right panel of Figure [2] shows the 
product of the reconstruction and identification efficiency (points) and the total efficiency (solid lines). 

Several sources of systematic uncertainty have been estimated by varying the assumptions applied to the analysis: 
the definition of tight and non-tight cuts, the event counting procedure, the photon energy scale, and the photon energy 
resolution. The uncertainties are dominated by the variations in the choice of the photon identification cuts and the 
precise choice of isolation cone properties, and a pj- and centrality-independent uncertainty of 31% is applied. 

1 dN N sis 

The photon yield is defined as j^jp^ (pi,c) = €[ - xN A xA p T where N A is the background-subtracted yield in 
region A, e tot is the abovementioned total efficiency, N evt is the number of events in the centrality bin c, Apj is the width 
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Figure 2: (left) The factor 1 — P (where P is the purity) extracted from data in each pj interval for four centrality bins, using the double 
sideband method |5|. (right) Photon reconstruction efficiency as a function of photon pj and event centrality averaged over \n\ < 1.3, based on 
MC calculations performed in four subintervals in r\. The points indicate the identification efficiency, without the contribution from the isolation 
requirement. The solid lines indicate the smoothed total efficiency used to correct the measured photon yield |'5'|. 



of the transverse momentum interval. The final yields are shown divided by the mean nuclear thickness (Taa) (which 
scales as the number of binary collisions) as a function of photon pj in Figure [3] While the most peripheral interval 
(40-80%) is unsealed, each subsequent centrality interval is scaled up by a factor of 10 for visibility. Also included are 
CMS data [4] for 0-10% central heavy ion collisions, \rj\ < 1.44 (a 10% larger interval in //), and an isolation condition 
of at most 5 GeV transverse energy in a cone of radius R[ so = 0.4, superimposed on the most central bin, as well as 
CMS pp data |4| from -\fs - 2.76 TeV, superimposed on the most peripheral bin. The systematic uncertainties on the 
yield are indicated by yellow bands, and the uncertainties on the mean nuclear thickness are tabulated in Ref. [5]. The 
ratio of the data to JETPHOX 1.3.0 ITBI [171 is shown in the right panel of Figure [3] It is observed that, within the 
statistical and systematic uncertainties, the data are consistent with the predictions in all centrality intervals. 



6. Conclusion 

Yields of isolated prompt photons in lead-lead collisions have been measured as a function of collision centrality 
in a kinematic range of photon pi — 45 - 200 GeV and for \q\ < 1.3. Photons have been reconstructed using 
the large-acceptance, longitudinally segmented ATLAS electromagnetic calorimeter, after subtraction of the average 
background, event-by-event. After scaling the yields by the mean nuclear thickness (Taa), they are observed to be 
constant as a function of centrality within experimental uncertainties, implying a linear scaling with the number of 
binary collisions. The scaled yields, as a function of pj, are found to be in good agreement with next-to-leading order 
pQCD calculations, as implemented in JETPHOX 1.3.0. 



References 

[1] ATLAS Collaboration, Phys. Rev. Lett. 105, 252303 (2010). 

[2] X. -N. Wang, Z. Huang and I. Sarcevic, Phys. Rev. Lett. 77, 231 (1996). 

[3] S. S. Adler et al. [PHENIX Collaboration], Phys. Rev. Lett. 94, 232301 (2005). 

[4] CMS Collaboration, Phys. Lett. B 710, 256 (2012). 

[5] ATLAS Collaboration, ATLAS-CONF-2012-051 (2012). 

[6] ATLAS Collaboration, JINST 3, S08003 (2008). 

[7] ATLAS Collaboration, [arXiv-T2 08.1967 [hep-ex]. 

[8] ATLAS Collaboration, |arXiv;090 1 ,0512| [hep-ex] . 

[9] https : //twiki . cern. ch/twiki/bin/view/AtlasPublic/LuminosityPublicResults#2011_Pb_Pb_Heavy_Ion_Collisions 

[10] T. Sjostrand, S. Mrenna and P. Z. Skands, JHEP 0605, 026 (2006)! 

[1 1] X.-N. Wang and M. Gyulassy, Phys. Rev. D44, 3501 (1991). 

[12] S. Agostinelli et al. [GEANT4 Collaboration], Nucl. Instrum. Meth. A 506, 250 (2003). 

[13] ATLAS Collaboration, Phys. Lett. B 707, 330 (2012). 

[14] ATLAS Collaboration, Phys. Rev. D 83, 052005 (201 1). 

[15] ATLAS Collaboration, Phys. Lett. B 706, 150 (201 1). 



/Nuclear Physics A 00 (2012) 745] 



5 



> 

a> -i n 7 



J3 
Q. 



h- 



Q. 
"D 



< 1Q 5 



no 3 



— 10 



10' 



10"' 



ATLAS Preliminary 
Pb+Pb^ = 2.76TeV 
V l_i nt = 133 lib" 1 

y hi<i.3 

y 



■ 

— i ■ 



- 0-10%(x10 ! ) 

- 10-20%(x10') 

- 20-40%(x10°) 

- 40-80% 
JETPHOX 1.3.0 

■ PYTHIAMC11 

- CMS HI (M<1.44) 

- CMS pp (h|<1 .44) 




100 



200 

photon p T [GeV] 



x 
o 

X 
CL 
H 
LU 

Q 



H^-f-^J _ 0-10%(+6) 

6 

i+ ^+ , JT 10-20%(+4) 

■ r + ; 

j,* I T ^ 20-40%(+2) 

2 

----±Jt± " 80% 

T TT y~ ATLAS Preliminary 

Pb+Pb ^ = 2.76 TeV - 
L. = 133^ib' (ll<1.3 

o 

100 200 300 

photon p T [GeV] 



Figure 3: (left) Efficiency corrected yields of prompt photons in \rj\ < 1.3 using HI tight cuts, isolation cone radius R; s0 = 0.3 and isolation energy 
of 6 GeV (5). Statistical errors are shown by the error bars. Systematic uncertainties on the photon yields are combined and shown by the yellow 
bands, (right) Efficiency corrected yields of prompt photons in < 1.3 using HI tight cuts, isolation cone radius i?j so = 0.3 and isolation energy 
of 6 GeV, divided by JETPHOX 1.3 predictions, which implement the same isolation selection. Statistical errors are shown by the error bars. 
Systematic uncertainties on the photon yields are combined and shown by the yellow bands. 
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